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Abstract — An adaptive joint power allocation (JPA) and linear 
receiver design algorithm using the minimum bit error rate 
(MBER) criterion for a cooperative Multiple-Input Multiple- 
Output (MIMO) network is proposed. The system employs 
multiple relays with Distributed Space-Time Coding (DSTC) 
schemes and an Amplify-and-Forward (AF) strategy. It is de- 
signed according to a joint constrained optimization algorithm 
to determine the MBER power allocation parameters and the 
receive filter parameters for each transmitted symbol. The 
simulation results indicate that the proposed algorithm obtains 
performance gains compared to the equal power allocation 
systems and the minimum mean square error (MMSE) designs. 

I. Introduction 

MIMO communication systems employ multiple collocated 
antennas at both the source node and the destination node in 
order to obtain the diversity gain and combat multi-path fading 
in wireless links. The different methods of STC schemes, 
which can provide a higher diversity gain and coding gain 
compared to an uncoded system, are also utilized in MIMO 
wireless systems for different sizes of the transmitter and dif- 
ferent conditions of the channel. Cooperative MIMO systems, 
which employ multiple relay nodes with antennas between the 
source node and the destination node as a distributed antenna 
array, apply distributed diversity gain and provide copies of 
the transmitted signals to improve the reliability in the wireless 
communication system |4|. Among the links between the relay 
nodes and the destination node, a cooperation strategy, such as 
Amplify-and-Forward (AF), Decode-and-Forward (DF), and 
Compress-and-Forward (CF) |l2] and various DSTC schemes 
in |3|, |22| and |23| can be employed. 

Since the benefits of the cooperative MIMO systems are 
noticed, extensive studies of the cooperative MIMO networks 
have been operated fSl-fT?!. In fTO'l, an adaptive joint relay 
selection and power allocation algorithm based on the MMSE 
criterion is designed. After each transmission, the optimal re- 
lay nodes with minimum MSE will be selected and the power 
allocation vectors will be fed back to the relay nodes. A joint 
transmit diversity optimization and relay selection algorithm 
for the DF cooperating strategy is designed in [ 11 1. A transmit 
diversity selection matrix is introduced at each relay node in 
order to achieve a better MSE performance by deactivating 
some relay nodes. In |8|, an optimal design for the amplify 
matrix at the relay node using the AF cooperating strategy is 
derived. The simplest three-node cooperative MIMO network. 



which contains one source node, one destination node and only 
one relay node in the system, is employed. The algorithm 
in |8| indicates that the optimal amplify matrix for AF 
strategy depends on the channel matrices when the direct Unk 
between the source node and the destination node is ignored. 
In |9|, a selection algorithm for one relay node among the 
other neighboring relay nodes is derived. However, with the 
simplicity of selecting only one relay node to construct the 
three-node cooperative MIMO network, the diversity gains of 
the system will be sacrificed. A central node which controls 
the transmission power for each link is employed in 1. 12] . 
Although the central control power allocation can improve the 
performance significantly, the complexity of the calculation 
increases with the size of the system. The works on the 
power allocation problem for the DF strategy measuring the 
outage probability in each relay node with a single antenna 
and determining the power for each link between the relay 
nodes and the destination node, have been reported in |15]-|EI- 
The diversity gain can be improved by using the relay nodes 
with multiple antennas. When the number of relay node is the 
same, the cooperative gain can be improved by using the DF 
strategy compared with a system employing the AF strategy. 
However, the interference at the destination will be increased 
if the relay nodes forward the incorrectly detected symbols in 
the DF strategy. 

In this paper, we propose a joint MBER adaptive power 
allocation algorithm with linear receiver design for cooperative 
MIMO systems. There are multiple relay nodes with multiple 
antennas to achieve an AF cooperating strategy with DSTC 
between the source node and the destination node. By using 
adaptive algorithms with the MBER criterion ||251 ||261 . the 
power allocation parameters and the linear receive filter pa- 
rameters can be determined and fed back to each transmission 
node through a feedback channel that is assumed error free 
and delay free. The performance indicates the advantages in 
diversity gain of the proposed JPA algorithm. 

The paper is organized as follows. Section II provides 
the two-hop cooperative MIMO system with multiple relays 
applying the AF strategy and DSTC scheme. Section III 
describes the constrained power allocation problem and the 
linear MBER detection method, and in Section IV, the pro- 
posed iterative SG algorithm is derived. Section V focuses 
on the results of the simulations and Section VI contains the 
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Fig. 1. Cooperative MIMO System Model with rir Relay nodes 



conclusion. 



II. Cooperative System Model 



The communication system under consideration, shown in 
Fig.l, is a MIMO communication system transmitting through 
a MIMO channel from the source node to the destination node. 
The BPSK modulation scheme is employed in our system by 
mapping the information bits b[i] to the transmitted symbols 
s[i] G {±1}- There are ti^ relay nodes with N antennas 
for transmitting and receiving, applying an AF cooperative 
strategy as well as a DSTC scheme, between the source node 
and the destination node. A two-hop communication system 
that transmits symbols from the source to rir relay nodes as 
well as to the destination node in the + 1 phases, followed 
by transmitting the amplified and re-encoded symbols from 
each relay node to the destination node in the next phases. 
We consider only one user at the source node in our system 
that has N Spatial Multiplexing (SM)-organized data symbols 
contained in each packet. The received symbols at the k — th 
relay node and the destination node are denoted as rsR,. 
and rsD, respectively, where k = 1,2,..., rir- The received 
symbols vsr^ will be amplified before mapped into an STC 
matrix. We assume that the synchronization in each node is 
perfect. The received symbols at the destination node and each 
relay node can be described as follows 



(1) 



=Ds\i]As\i]s\i] + ns\i], 

1^1,2, ... ,N, 

where the diag( ) operator transforms the argument into 
a diagonal structure and the (n^ + l)N x 1 vector ns[i] 
denotes the zero mean complex circular symmetric Additive 
White Gaussian Noise (AWGN) vector with variance a^. 
The transmitted symbol vector s[i] contains N parameters, 
s[i] — [si[i], S2[i], SN[i\]'^ , which have a covariance ma- 
trix i?[s[i]s-^[i]] — (j'^In, where E[-] stands for expected 
value, (•)^ denotes the Hermitian operator, cr^ is the signal 
power which we assume to be equal to 1 and In is the 
N X N identity matrix. The diagonal matrix -D5[i] is the 
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{rir + 1)^ X rirN channel gain matrix that contains the 
N X N channel matrix H[i] between the source node and 
the destination node and the channel matrix Fk[i] between 
the source node and each relay node, where k = 1,2, ...,nr. 
The rirN x N matrix As[i] denotes the power allocation 
matrix for the transmitted symbol vector s[i]. Notice that the 
{nr + l)N xl vector rs[i] in ([T]i contains the received symbol 
vector at the destination node and the relay nodes, which is 
rs[i\ = [rsD[i\,rsRAi\,---,rsR„^[i]]'^. 

After processing the received vector rsR^[i] at the k — th 
relay node, the signal vector Ss^^. [i] can be obtained and 
will be forwarded to the destination node. A power allocation 
parameter vector will be assigned to ssr^ [i] first, and then the 
amplified symbols in ssr^ [i] will be re-encoded by a x T 
DSTC scheme M{s) and then forwarded to the destination 
node. The relationship between the k — th relay and the 
destination node can be described as 



k^l,2. 



(2) 



where the TV x T matrix Mr^dIi] is the DSTC matrix 
employed in the relays whose elements are the amplified 
symbols in ssR^i]. The N x T received symbol matrix 
RRkoii] in © can be transformed and denoted as a NT x 1 
vector rji^jjli] given by 



rR,D[i\ = G',,[i\AR^D[i\ssRAi\ + '>^RkD[i\ 



(3) 



1,2, 



X N matrix An^oii] = 
contains the power allocation 



where the diagonal N 
diag[ai^^^[i],...,ajvH^„[z]] 
parameters assigned for the j — th symbol The 
NT X N matrix G'j. [i] stands for the equivalent channel matrix 
which is the DSTC scheme M{s[i]) combined with the 
channel matrix GR^oii]- The NT x 1 noise vector nji^uli] 
is an equivalent noise vector at the destination node which 
contains the noise parameters in Np^^jjli]. After rewriting 
-RflfccH we can consider the received symbol vector at the 
destination node as n,. + 1 parts, one is from the source node 
and the remaining rir are from the relay nodes, and write the 
received symbol vector for cooperative detection as 



H[i]AsD[i]s[i] 
G[\i]Aji^D[i]ssRi [i\ 



r I = 



G'2\i]AR^D\i]ssR2 [i] 
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(4) 



Dd [i] Ad [i]sD [i] + no [i] , 



where the (7irr+l)Afx (n^+ljiV diagonal matrix JDo [z] con- 
tains the channel gain of all the links between the relays and 
the destination node. The {rir + 1)^ x N power allocation pa- 
rameter matrix AD[i\ = [AsD[i\,AR^D[i\, ■■■ , Ar^^dIiW 
is constructed by all the power allocation matrices between all 
the relay nodes and the destination node. 



III. Joint linear MBER receiver design with Power 
Allocation 

The MBER receiver design with power allocation for 
every link in all the phases between the transmitting 
nodes and the receiving nodes is derived as follows. If 
we define a {rirT + 1)N x 1 parameter vector Wj[i] ~ 
[wj-^[i],Wj2[i],...,Wj^^^-^[i]] for i — 1,2,...,N to determine 
the j — th symbol Sj[i], by using (4), the desired information 
symbols at the destination node can be computed as 



sgn(u)f [i]r[i]) = sgn(sj[i]), 



(5) 



where sgn(-) is the sign function and Sj[i] denotes the detected 
symbol at the receiver which can be further written as 



Si I 



wf [i]r[i] 



= wf\i]{DD\i]AD\i]sD[i\ + noii]) 

„ „ (6) 

= wf\i]DD[i\AD\i]sD[i\+wf[i\nD\i] 

where s'j[i] is the noise-free detected symbol, and ej[i] 
denotes the error factor for the j — th detected symbol s'j [i] 
with zero mean and variance a'^w^[i]wj[i], where cr^ is the 
variance of the received noise Define a N x Nj, matrix 

S which is constructed by a set of vectors s and Ni, = 2^, 
containing all the possible combinations of the transmitted 
symbol vector s[i] and we can obtain 



Sj\i] = wf[i\r[i\ 



wf[i]DD[i]AD[i]sS] 
; = l,2,...,iVfc, 



(7) 



where Sj [i] is an element from the noise-free detected vector 
si[i]. To derive the BER expression for the linear filter Wj[i] 
and the power allocation parameter matrix we should 

determine the probability density function (pdf) of f[i] which 
is 



1 



Pr[i] =- 



E 

1=1 



exp( 



(Sj[i]~sgn(bj[i])sj[i])^ 
2a>f[i]tOj[i] ' 



(8) 



As a result we can obtain the bit error rate (BER) expression 
as 



Pe {wj \i] , aj [i] ) = — ^ Q (cj \i] ) 



1=1 



where 



sgn{sj\i])sj[i\ 
an\ w^\i]w,\i 



(9) 



(10) 



Before we derive the adaptive algorithm, we have to spec- 
ify the power allocation parameter aj[i] in (|9]l. The power 
allocation matrices, A5£i[i], Ag/j^, [i] and A]^^£)[i], where 
k — l,2,...,nr, have to be estimated in our algorithm. For 
simplicity of expression, we use the factor a in our estimation 
function instead of the power allocation matrices. It is shown 
that in what follows, the parameters a 



]SD I 



and 



'^3r^d[A denote the power allocation parameters assigned for 
the j—th signal in A.s_D[i], Asr^. [i] and Afl^^H, as described 
by 



AsD [i\ = diag . , Asr^ [i] diag 



AB.^D[i] = diag 



By substituting ([8]l and ( fTOl ) into (|9]l and taking the gradient 
with respect to different arguments, we can obtain (flU - (fT4b . 
where hj[i] is the j — th column of H[i\ and QnkD-lA ™d 
f SR.^ . [*] denote the j — th column of Gi^^^/jH and FsRk 
respectively. By making use of (fTTT i - (fl4] i and an SG algorithm 
in f\\, the updated receive filter Wj [i] and the power allocation 
parameters assigned for the j — th transmitted symbolsj[j] 
between the source node and the destination node aj^^j [i], 
between the source node and the k — th relay node cnj^j^^ [i], 
and between the k — th relay node and the destination node 
'^jn r>M' '^aii obtained as follows 
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(15) 


ajSD [« + !] = 






(16) 


"jr,d[« + 1] = 






(17) 


"isfl, [« + !] = 






(18) 



where /i and 7 are step sizes for estimating the receive filter 
parameters and the power allocation parameters, respectively. 

A. Adaptive MBER SG Estimation and Power Allocation 

The key to achieving the adaptive SG estimation algorithm 
described in (fTsTl-lfrsIl is to find out an efficient and reliable 
method to calculate the pdf of the received symbol vector 
r[i] at the destination node. In I30l - li32i . the Parzen window 
and kernel density estimation method are introduced, which 
can guarantee the accuracy of the probability distribution. 
Because the noise vector in (4) is Gaussian and the kernel 
density estimation can provide the reliable pdf estimation 
when dealing with Gaussian mixtures, we choose the Kernel 
density estimation method in our study. 

By transmitting a block of K training samples s[i] — 
sgn(b), the kernel density estimated pdf of s[i] is given by 



Ps[. 



KV27rpnJwf\i]w^ 



3 

( _ jsjli] - sgn{bj\i])sj\i])^ . 



(19) 



where p„ is related to the standard deviation of noise (t„ and 
it is suggested in |30| that a lower bound of /?„ = (3^) ^ On 



1 c2[jl sgn{sjl{])rli]a„Jwf[{\wj[i] - sgn{sj[{\)wf [{\r[{]an{wf [{]wj[i]) lwj[i] 
S/Pe^ [i\ = ^Yexp( - 



(11) 



(12) 



1 ^ c2[i] s3n(sj[i])'ii)|^gfl^o,.Haji, bH/sh^.W^jW 
vPb„ j = — = > ex'p ^ — ) — 



(13) 



(14) 



should be chosen. The expression of the BER can be derived 
as 



where 



1=1 

sgn{sj[i\)sj[i] 



Pr 



[i\wj [i] 



(20) 



(21) 



Then by taking the gradient of (20) with respect to different 
arguments we can obtain the gradient value of the estimated 
error probability, and by iterative calculation of the parameters 
in the receive filter vector with the power allocation parame- 
ters, we can finally achieve the optimal joint estimation result. 

B. Adaptive SG Channel Estimation 

In this subsection we will derive an adaptive SG algorithm 
for estimating the equivalent channel matrix Do[i] - The chan- 
nel estimation can be described as an optimization problem 

Dd[i] - arg min E[\\r[i] ~ D D[i]A[i]DSDm'']- (22) 

Dd[i\ 

Define the received symbol vector r[i] = [ri[i], r„,,-)_i [i]]. 
The optimization problem described in (l22t can be divided into 
three parts, which correspond to individually computing H[i\ 
and G'j^^d[i\. We then derive an SG algorithm by calculating 
the cost function Gh and Cg'^ 

CH = E[\\r^[i]-H\i]AsD[i]s\i]\W (23) 
Cg', ^E[\\ru\A^ G'^\i]AR,Dmm^]- (24) 

Then by taking instantaneous gradient terms of Gh and Cg^ 
with respect to and the j — th column of the equivalent 
channel matrix G/j^ we can obtain 

VCh.[,] = -AsD[i]s\iKri\i] - H[^]AsD\^]s[^]), (25) 

^Cc'im =-Afc^[z]s[i](rfcH-Gij^^[z]Afl,,z3[z]s[z]). (26) 

Considering the SG descent rules O] and the results of gradient 
terms we can obtain the adaptive SG channel estimation 
expressions which are 

H[t + l]=H\i]~pVCH'ii] (27) 



G'k[i + l]^Gi[i]~/3VC, 



G'tl 



(28) 



where /3 is the step size of the recursion. The adaptive SG 
algorithm for the equivalent channel H[i] and Gjj[i] requires 
the calculation complexity of {0{N^)) and (O(TiV)). 

IV. Simulations 

The simulation results are provided in this section to assess 
the proposed algorithm. The BER performance of the adaptive 
JPA and the EPA using a linear MBER receive filter (JPA- 
MBER) and (EPA-MBER) algorithms with the power con- 
straint employs different numbers of relay nodes and different 
DSTC schemes in f24| are compared with — 1,2 relay 
nodes in the simulation. In the simulation we define the power 
constraint Pt as equal to 1 , and the number of antennas N ~ 2 
at each node. In Fig. 2, the Almouti STBC scheme is used at 
relay nodes in the simulation. The results illustrate that the 
performance of the EPA-MBER and JPA-MBER algorithms 
is superior to that of MMSE, and the JPA-MBER algorithm 
outperforms the EPA-MBER algorithm using the same DSTC 
scheme to achieve the identical BER. There is about 2.5 dB 
of gain compared the JPA-MBER algorithm with the EPA- 
MBER one when using 1 relay node to achieve the BER in 
10^'^. The performance improvement of the proposed JPA- 
MBER algorithm is achieved with more relays employed in 
the system as an increased spatial diversity is provided by 
the relays. The coding gain achieved by using different DSTC 
encoders at the relay nodes, namely, the Alamouti and the 
randomized Alamouti (R-Alamouti) is shown in the results 
of Fig. 3. The number of relay nodes is equal to 2. The 
JPA-MBER algorithm outperforms the EPA-MBER algorithm 
when using a different STC scheme, and the JPA-MBER 
algorithm with R-Alamouti STBC scheme obtains 1 dB of 
gain compared to the that of the system with the Alamouti 
STBC scheme. The proposed adaptive algorithm can be used 
in the cooperative MIMO system with linear detector in ifTsl 
or non-linear detector in [17 1. With different criteria in [l27ll 
to f29l, the proposed adaptive JPA algorithm can obtain more 
diversity gains compared with the EPA algorithms. 

V. Conclusion 

We have proposed a joint power allocation and receiver 
design algorithm using a linear MBER receive filter with the 
power constraint between the source node and the relay nodes, 
and between relay nodes and the destination node. A joint iter- 
ative estimation algorithm for computing the power allocation 



parameters and linear MBER receive filter vector has been 
derived. The simulation results illustrate the advantage of the 
proposed power allocation algorithm by comparing it with the 
equal power algorithm. The proposed algorithm can be utilized 
with different DSTC schemes and a variety of detectors [17] 
|[T8| and estimation algorithms lIZTl ||28| in the cooperative 
MIMO systems with AF strategy and can also be extended to 
the DF cooperation protocols. 
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